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Abstract 

The key challenges for models with large extra dimensions, posed by neutrino 
physics are: first to understand why neutrino masses are small and second, 
whether one can have a simultaneous explanation of all observed oscillation 
phenomena. There exist models that answer the first challenge by using sin- 
glet bulk neutrinos coupled to the standard model in the brane. Our goal in 
this paper is to see to what extent the simplest versions of these models can 
answer the second challenge. Our conclusion is that the minimal framework 
that has no new physics beyond the above simple picture cannot simultane- 
ously explain solar, atmospheric and LSND data, whereas there are several 
ways that it can accommodate the first two. This would suggest that confir- 
mation of LSND data would indicate the existence of new physics either in 
the brane or in extra dimensions or both, if indeed it turns out that there are 
large extra dimensions. 

PACS: 14.60.Pq; 14.60.St; ll,10.Kk; 



I. INTRODUCTION 

Particle physics models where there are large hidden space dimensions beyond the three 
familiar ones have been the focus of intense activity during the past two years . Beyond 
the simple reason that such extra dimensions are predicted by string theories, a major point 
of interest in these models is that often these large extra dimensions come with a TeV scale 
for the strings which leads to a plethora of new observable phenomena in collider as well in 
other arenas of particle physics and cosmology. 
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In order for these models to provide a satisfactory description of low energy particle 
physics, they must handle some obvious problems that come with the existence of a funda- 
mental scale in the multi-TeV range. One such problem has to do with understanding the 
small neutrino masses in a natural manner. The conventional seesaw explanation which 
is believed to provide the most satisfactory way to understand this, requires that the new 
physics scale (or the scale of SU{2)f( x U{1)b-l) be around 10^^ GeV or higher. Clearly, low 
string scale theories do not have any fundamental scale of that type. Moreover, the low value 
of the string scale leads to enhanced (and unacceptable) contributions to neutrino masses 
from higher dimensional operators. While there are suggestions involving thick branes with 
point splitting ||^ to remedy similar problems that arise from AB ^ operators, they don't 
work for neutrino mass operators. Therefore, a necessary ingredient to understand small 
neutrino masses in the low string scale models is to assume that theory have a B — L sym- 
metry. This will forbid higher dimensional operators LHLH/M* {M* is the string scale), 
which are the source of the problem for neutrino masses. 

Depending on whether i? — L is a global or local symmetry, one can have two ways to 
solve the neutrino mass problem in models with large extra dimension^. In the former case, 
discussed in 0, one has to introduce singlet bulk neutrinos which then lead to small Dirac 
masses for them. On the other hand, if B — L symmetry is chosen to be a local symmetry, 
anomaly cancellation requires that, right handed neutrinos be present in the brane as in the 
models discussed in Ref. Since local B — L must be broken to avoid massless gauge 

bosons, one again has to deal with the induced operators of the same type as above with 
M* = Mb^l- It was shown in Ref. |]^,|| that to get neutrino masses in the desired eV range, 
one must have string scale M* ~ Mb-l ~ 10^ GeV range or higher. As a result, a class of 
experimentally accessible phenomenological predictions are lost, although long range gravity 
tests are still possible. These models are similar to the ones discussed in 

In the context of models that have global U{1)b-l symmetry, one can maintain the TeV 
scale for the strings and still get small neutrino masses by introducing isosinglet neutrinos 
in the bulk as has already been discussed in Ref. . These models are interesting because a 
very minimal set of particles beyond the standard model are sufficient to get small neutrino 
masses. Key reason for this result is the relation between the fundamental scale, M*, the 
radius of the extra dimensions, i?, and effective Planck scale, Mp^, 

Since the singlet neutrino is a bulk field, the effective couplings of its Fourier modes to the 
standard model fields are naturally suppressed by the ratio M* /Mpe, which for a TeV M* 
produces the right order of magnitude for neutrino masses. 

An analysis of the implications of the mixing profile in these models for solar neutrino 
deficit was discussed in |T^. Also, implications for atmospheric neutrinos were discussed 
in [|lT|, and some phenomenological bounds were given in |ToHT3|. Our goal in this paper is to 



^It is generally believed that string theories do not have any global symmetries which would 
seem to imply that B—L must be a local symmetry. In our discussion, however, we will consider that 
B — Lis a. global symmetry of the theory, as in Ref. and see where it leads us phenomenologically. 
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attempt a unified explanation of all known neutrino oscillation data i.e. solar, atmospheric 
as well as LSND under different assumptions for the initial input parameters for the minimal 
bulk neutrino scenario. 

The basic reason for embarking on such an ambitious program is the encouraging feature 
that the masses of the lower KK modes of the bulk neutrinos are given by integral multiples 
of which is of order 10~^ eV when R ~ millimeter. This is of the right order necessary to 
solve the solar neutrino problem via small angle MSW mechanism using z/g to oscillation. 
Thus we see that there exists a natural way to understand the lightness of the sterile neutrino 
Pl jnilJT^ , a situation if realized would pose a major challenge to four dimensional theories. 
Once the solar neutrino problem is understood, it would appear that all the necessary 
ingredients are at hand to understand the atmospheric and LSND data using oscillations 
among the familiar neutrinos i.e. z/g for LSND and for atmospheric. 

This program has already been undertaken in special parameter domains and it has 
already been suspected [|,|n|] that it does not really work. Our goal is to extend these 
discussions to a somewhat larger parameter domain to see if there is chance for this program 
to succeed and unfortunately our answer is also in the negative. Our work complements the 
above works and extends them. Specifically, we try to give analytical reasonings to see how 
the different oscillation data can (or cannot) be understood. Since we do not take recourse to 
a detailed numerical analysis, we cannot rule out the possibility that some small parameter 
domain exists where all data can be accommodated; but we consider that unlikely. 

The negative conclusion of our work, combined with the works of implies with 

virtual certainty that in the large extra dimension framework, understanding the neutrino 
data does require new physics beyond the standard model in the brane or new physics in 
extra dimensions or both. 

Our basic strategy is as follows: we start by requiring that the parameters of the theory 
provide an explanation for solar and atmospheric oscillations. Then we ask whether they 
can account for the small LSND probability of z/g appearance from the z/^ beam. 

This paper is organized as follows: in section II, we discuss neutrino oscillations with 
a single flavour, which depends on two parameters, the bulk radius and a dimensionless 
parameter ^ related to the Dirac mass term in the theory. The latter defines the pattern of 
neutrino oscillations. This analysis sets the stage for the three flavour case, which we discuss 
in section III for various possible domains of the parameter space. We end the paper with 
a concluding section that summarizes the results and discusses their implications. 



II. OSCILLATIONS WITH A SINGLE FLAVOUR 

Let us begin our discussion by focusing on the simplest case with one generation of 
fermions in the brane and one bulk neutrino, to understand the general profile of the neutrino 
oscillations in models with large extra dimensions. We will discuss the necessary ingredients 
to understand the three flavour case that is explored in the next section. Obviously, the 
fields that could propagate in the extra dimensions are chosen to be gauge singlets. Let 
us denote bulk neutrino by i>B{x^,y). It has a five dimensional kinetic energy term and a 
coupling to the brane field L{x^) given by 

C = kLHubr{x, y = 0)+ dy ubl{x, y)d5UBR{x, y) + h.c. (2) 
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where from the five dimensional kinetic energy, we have only kept the 5th component that 
contributes to the mass terms of the KK modes in the brane; H denotes the Higgs doublet, 
and K = h^^^ the suppressed Yukawa coupling. It is worth pointing out that this suppression 
is independent of the number and radius hierarchy of the extra dimensions, provided that 
our bulk neutrino propagates in the whole bulk. For simplicity, we will assume that there 
is only one extra dimension with radius of compactification as large as a millimiter, and the 
rest with much smaller compactification radii. The smaller dimensions will only contribute 
to the relationship (|l]) but its KK excitations will be very heavy and decouple from neutrino 
spectrum. Thus, all the analysis could be done as in five dimensions. 

The first term in Eq. (|^) will be responsible for the neutrino mass once the Higgs field 
develops its vacuum. The induced Dirac mass parameter will be given by m = /tf , which 
for M* = 1 TeV is about h ■ 10~^ eV. Obviously this value depends only linearly on the 
fundamental scale. Larger values for M* will increase m proportionally. After introducing 
the expansion of the bulk field in terms of the KK modes, the Dirac mass terms in could 
be written as 

/- -/ ^ ( rri ^/2m \ ( uqb \ ro\ 
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where our notation is as follows: v'^ represents the KK excitations, the off diagonal term 
^J2m is actually an infinite row vector of the form A/2m(l, 1, ■ ■ ■). The operator stands 
for the diagonal and infinite KK mass matrix whose n-th entrance is given by njR. This 
notation was introduced in [|^ to represent the infinite mass matrix in a compact manner. 

Using this short hand notation makes it easier to calculate the exact eigenvalues and the 
eigenstates of this mass matrix Simple algebra yields the characteristic equation 

2A„ = 7re'cot(7rA„), (4) 



with A„ = m„i?, ^ = \/2mR, and where m„ is the mass eigenvalue |]6|,p!0[|. The eigenstates, 
on the other hand are given symbolically by p| 



1 

' r). 



(5) 



where the sum over the KK modes in the last term is implicit. is the normalization 
factor given by 

Nl = ^ (A^ + /(O) , (6) 

where /(^) = ^^/2 + 7r^^'^/4. Figures 1 and 2 depict the exact numerical results of A„ and 
Nn for several choices for the value of the ^ parameter. Using the expression (^), we can 
write down the weak eigenstate in terms of the massive modes as 

OO 

n=0 

Thus, the weak eigenstate is actually a coherent superposition of an infinite number of 
massive modes. Therefore, even for this single fiavour case, the time evolution of the mass 
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eigenstates involves in principle all mass eigenstates and is very different from the simple 
oscillatory behaviour familiar from the conventional two or three neutrino case. The time 
dependent survival probability is given by 



PsurviL) = \{uL{L)\uLmt 
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It is clear that the survival probability depends strongly on the parameter ^, reflecting the 
universal coupling of all the KK components of vb with i>l in (H). Figures 3 and 4 show the 
profile of the survival probability obtained from the numerical solutions for three different 
values of ^ and is clearly very different from simple familiar oscillatory behaviour. However, 
to better understand these results, we will follow an analytical approach in what follows. 

It is simpler to consider the two limiting cases. First let us assume that ^ -C 1. As already 
known PJTOd, the eigenvalues in this case are given by Aq = mR, and \n = n otherwise. 
Therefore, to a good approximation, we may take the mixing parameters A^o = 
Nn = {n/^)ri for non zero n. Where the extra factor 77 = (1 + vr^^^/G)^/^ is introduced to 
keep the proper normalization in the expansion (|^. This approximation is confirmed by our 
figures 1 and 2. The survival probability is now given as 
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It is simple to see from last expression that the probability has an oscillation length Lose = 
AnER^. A typical profile in this case is depicted in figure 4. Also, the main contribution 
to the oscillation pattern comes from the lowest elements of the tower, which turns out to 
be the main component of vl- Another way to see this result is to note that ^ could be 
made small by making R small; but this makes the KK excitation masses large so that they 
decouple from the light sector of the theory leaving only the right handed component of 
the zero mode of the bulk neutrino, vbr to form a Dirac mass with and stay light. The 
left-handed zero mode decouples and remains massless. 

In the <C 1 limit, the small mixing with the tower elements leads to an oscillation 
pattern dominated by the lightest KK mode with a mass just about 1/R and more or less 
simulate the familiar oscillatory one. This can happen, for instance, if ~ 0.2 mm, which 
gives ~ 10~^ e^^, just about what it is needed to provide an explanation to the solar 
neutrino problem assuming a small MSW mixing angle []TU|. To use this case to compare 
with solar neutrino data, one needs to include the matter effect, which has already been 
discussed in Ref. and we do not enter into this here. 

A more direct application of the formula in Eq. (9) can be made to discuss the atmo- 
spheric neutrino oscillations, which is a vacuum oscillation. However, to get the equivalent 
of large mixing angle we need to adjust ^ ~ 1 and l/i?^ ~ Am^^^, which means that Eq. 
(9) is not applicable and we must therefore use Eq. (8) and study its large ^ limit. 

To discuss this, we start with the limit when 1^ ^ 1. First thing to note is that the 
pattern of eigenvalues is very different in this case from the small case. For small values 
of n the eigenvalues are well approximated by A„ = while A^„ becomes n independent 
until certain cut off value that roughly speaking is given by tia = 7r^,^^/4. Beyond that 
point A„ ~ and we may no longer neglect the contribution of A„ to (Bl). This results in 
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a suppression of Nn that goes like 1/n. It should therefore be reasonable to consider as a 
first approximation that the expansion (|^) is cut off at and that all the mass eigenstates 
contributing to are equally suppressed by l/^/nJi. In this approximation, we obtain the 
survival probability to be 
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Certainly, this relation presents a oscillatory profile, however, the superposition of the equally 
suppressed oscillations will result most of the time on a destructive interference, with the 
exception of the very sharp resonances that appear each time L reaches a multiple of the 
oscillation length. In other words, is a superposition of a large number of mass eigenstates, 
with masses covering a large range, from 1/2 to r^A in units of 1/-R, all of them contributing by 
the same amount. As a result, once ul is released, one may surmise that the time evolution 
of the different components will most likely wash out the original coherent superposition 
and the initial ul will almost disappear. As fig. 3 shows, this conclusion is borne out by the 
the numerical analysis. To get an analytical result that also supports this conclusion, note 
that in the sum in Eq. (10), the dominant contributions come when n k in which case 
each term in the sum averages to 1/2 and on performing the double sum it is easy to see 
that one arrives at Pgurv ~ l/'^A- 

If TT-A is very large, it suppresses the survival probability too much and cannot help in 
the understanding of the atmospheric data. So, clearly, if we wanted an understanding of 
the atmospheric data, we must assume smaller ^, perhaps values closer to one. In this case, 
truncating the sum in the expression for the survival probability in Eq. (^, cannot be 
justified and we must seek an alternative way to deal with Eq. (||). An approach suggested 
m in is to use a continuous approximation to the sum in Eq. (8), which leads to 
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where z = L/2EB?. Clearly, for large ^; f{S,) ~ 7r^,^^/4, and last expression in Eq. (11) 
simplifies. In order to study the dependence of Psurvi.^) on ^, we plot Pgurv in Fig. 3 in 
the continuous and discrete approximations. One thing that emerges is the non-oscillatory 
nature of the function as we exceed ^ = 1. We also see that the continuous hmit gives a good 
approximation for the slope (see figure 3), even for cases where ^ < 1 [ll|]. This expression. 



however, does not represent the survival probability in the limit ^ <^ 1, since in this case 
(i.e. ^ ^ 1), the Psurv has an oscillatory behaviour unlike the last term in Eq. (11) (see 
figure 4). In the overlap region close to ^ ~ 1, one may evaluate the second term in Eq. (11) 
in a slightly different way which leads to an expression for the survival probability as 



,{0 = p\Q + ^\0. 



(12) 



where we have introduced the new variable Q = 

tions p(C) = 1 - C(C) - SiC); and ^(C) = C(C) 
Fresnel integrals: 



zf{^) = {7r^y2R)^L/2E and the func- 
S{() with S and C the sine and cosine 
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We draw attention to the fact that ( involves not only the model parameters ^ and R but 
also the experimental variable z. Thus for a given model, the variable ( has different values 
for different oscillation experiments; for instance, for the atmospheric neutrino case, a typical 
value for C ~ if 1.2 X 10^ eV'^ 

The expression (|1^) is valid just before PsurviC) reaches the average, Psurv = 4/7r^^^. In 
figure 5 we present the behaviour of p, ^ and Psurv Note the steep fall off of Psurv as ( in- 
creases. For C = 1; Psurv has already dropped under 20%, which is smaller than the observed 
deficit in solar and atmospheric data. Therefore, if we want to fit the overall suppression of 
the atmospheric neutrinos, we must remain in a very narrow range of parameters. A rough 
estimate of these parameters may be obtained by expanding (|I2]) to first order on (, which 
yields 





Ps.™(C)^l-2^-C = l-20F(^) ( — ) . (13) 
Now, we may invert the last equation to get 



e ^ (1 - Pe.p) (±\ 

R JI^ \AeJ 
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(14) 



Using average values for the L,E and Pexp for atmospheric neutrinos, we find S,^/R ~ 
IQ-'^eV. Note that this naive approximation gives just about the value obtained by the 



numerical fitting of the data in |jTT]]. For this value, as ^ > 1, we get R in the range 
10"'^ eV < 1/R < 10^^ eV implying 1 < < 10. Thus a solution to the atmospheric 
neutrino puzzle using KK modes requires that R be in the sub-millimeter range. This in 
turn determines how large M* should be to avoid any unwanted large fine tuning. We see 
that for m > 10~^ eV and Yukawa coupling h of order one, we need M* > 100 TeV. 

We also point out that in the case of ^ > 1, extra dimensions must be much larger than a 
millimeter if we want to fit the solar neutrino data either via MSW or via vacuum oscillation. 
To see this note that if we take Am^^; ~ 10~^ eV^, this would imply L/4E ~ l/Am^^; ~ 10^ 
eV-2. Putting this in Eq. (13), we get ^'^/R ~ 10"^ eV. For ^ > 1, this implies that R > 0.2 
mm. Similar estimate for the case of vacuum oscillation yields R^ 2 cm. Therefore, in all 
our discussion of solar neutrino oscillations involving the bulk neutrinos, we will work in the 
approximation ^ <^ 1. 

Let us now summarize our findings for a single neutrino case: for ^ -C 1, Psurv has an 
oscillatory behaviour as given by (H). As ^ approaches 1, A„ starts to deviate from the 
integer value n, which in turn disturbs the periodic nature of Psurv and the maximum and 
minimum values of probability are not reached away from the source 0. This picture gets 
worse as one approaches ^ ^ 1 when a very sharp slope drives the probability near zero as 
we move away from the source of the neutrinos and it remains around its average, 4/7r^^^, 
most of the time. For this reason, whenever we try to fit the atmospheric neutrino data 
with bulk neutrinos, the value of ^ must be tuned to a very narrow range. 

Let us apply the discussions of this section to study how the oscillation of known neutrinos 
to bulk ones would effect the current experiments such as CHOOZ-PALOVERDE, LSND and 
the atmospheric neutrinos. For this purpose, we first note that in contrast with the usual two 
neutrino oscillation case, where the pattern is determined by three parameters, the mixing 
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angle 9, the Am^ and L/4E for the experiment, in the case of bulk neutrino oscillation, 
we have L/AE characterizing an experiment like before, bulk radius R, which replaces Am^ 
(which we will assume to be in the milli-meter range) and model parameter ^ (which is 
the analog of the mixing angle). For a given R, ^ is the only parameter characterizing the 
oscillation pattern. In Fig. 6, we plot the variation of the survival probability Psurv against ^ 
for the various cases mentioned above for a typical characteristic value of L/4E. We see from 
this figure that to explain the observed overall deficit of atmospheric neutrinos by nearly 
50%, one needs to go to ^ ~ 1.5 or so. This conclusion is in acoord with our conclusion 
based on Eq. (14) above. 



III. THREE FLAVOUR OSCILLATIONS 

Let us now apply the discussion of the previous section to the case of three standard 
model generations in the brane so that we have three brane neutrinos i^e,ti,T- To give masses 
to all of them in a minimal scenario, we will use three bulk neutrinos and allow arbitrary 
Yukawa couplings between the bulk and the brane neutrinos. This leads to an arbitrary 
Dirac mass matrix that involves the familiar left-handed neutrinos and the right handed 
Kaluza-Klein modes of the bulk neutrinos. As already discussed in the most general 
Dirac mass terms with three flavours may be written, after a rotation of the bulk fields, as 

C = vl-UMd ■ VBR{y = 0) + j dy vbl ■ ■ vbr + h.c, (15) 

where [/ is a unitary matrix and M/j = Diag{mi,m2,m3), in the basis where = 
ii^e,J^fi,^T)L', and ub = (z/]j, i/^, z/^). The mass parameters are just the eigenvalues of 
the Yukawa coupling matrix multiplied by v, the vacuum expectation value of the standard 
model doublet field. They are of the order of eV or less since the couplings of the bulk modes 
to the brane fields are naturally suppressed, as already stated on the previous section. 

Now, to simplify the discussion, we rotate the weak eigenstate neutrinos Ua to the ones 
related to the weak eigenstates by the rotation U i.e. z/^ = Uaa^a, where a = e,fi,T and 
a = 1, 2, 3. We then have 

3 

^ = E 

a=l 

This reduces the problem to a consideration of three KK towers mixing with three neutrinos 
which are related to the weak interaction eigenstates by the unitary matrix U defined above. 
Each tower is characterized by its parameter defined in analogy with section II as C,a = 
^J^rrtaR. After diagonalization of the mass matrices in each tower, each standard neutrino 
can be written as a coherent superposition of the three different towers of mass eigenstates: 

3 3 oo 

^a = Y Uaal^a = J2 UaaJ^i>ak- (17) 
a=l a=l k=0 

This expression generalizes Eq. (|^. It is now clear that the three flavour oscillations will 
correspond to the oscillations among the three towers. In this regards, the explanation to 
neutrino puzzles is not any more described in terms of three single neutrinos, as in the usual 



.L^BR^y = 0) + / 



dyVBLd^^BR 



+ h.c. 



(16) 
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case; instead all the KK modes can contribute (unless all the KK excitations decouple from 
the spectrum). 

To proceed further, let us define the partial transition probabilities 



Pap{L) ^ {ML)\u^{0)){up{L)\up{0)). (18) 

Notice that the diagonal component p^a may be interpreted as the survival probability of 
z/q,, and it takes the form of Eqs. @ and (0) in that case. It of course involves the u^s and 
not the flavor eigenstates as is obvious. The transition probability among standard flavours 
can be written in terms of the p^/j's as 

Pab = E U:^U,^U:^Uap Po^p. (19) 

Neglecting all CP phases we may expand Eq. (|T^) to the form 

sin^ ff (AL - AL)) 

Now, we are ready to address the oscillation problem. Our approach will be as follows: 
we will first select the parameter range that provides overall reduction required to explain 
the solar and atmospheric data, and then ask whether for the same range of parameters we 
can explain the observed oscillation between z/^ to z/g reported by LSND. Without loss of 
generality, we can assume the hierarchy < ^2 < ^3- We consider three possible scenarios: 

(i) 6,2,3 < 1; 

(ii) 6,2 < 1 < ^3 and 

(iii) 6 < 1 < 6,3- 

The case where all > 1 is already ruled out since, as we discussed earlier, it cannot 
explain the solar neutrino data without implying that the extra dimensions be too large. We 
therefore do not discuss this case. In the discussion of our results, the following expressions 
for the partial transition probabilities will be very useful: If S,a,i3 ^ 1, then 



Pa(3 — -5-^ 



COS (^Am^,) + {a + if) dz) + iie, {c\z) + s\z))\ , (21) 



with Am^^ = — m^; r^^ = 1 + and z = L/2EB? as before, and where we have 

introduced the functions 

( c(2;) \ ]_( cos(zn2) \ 



It is simple to check that the diagonal term of Eq. (|2lD reduces to Eq. (H). For <^ 1 < ^/j 

we get 

1 ^2 

Pa/3 = —Cp{z) + -| {c{z) Cp{z) + S{z) Sp{z)) ' (23) 
Va Va 
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cos(C|)p(C/3) - sin(C|)'i9(C,3) and spiz) = sm{Q)p{Cp) + 



where we now used C[^{z) — uuni^^^jy^i^i,/:); — oinvs/j; 

cos(C|)'*9(C/3) with (/^ = yi7r^^/2, and p and 'd as defined in the previous section, we should 
stress that all those functions have a similar deep behaviour. And, finally, for ,g > 1 we 
have 



Pa/3 = Ca{z) Cp{z) + Sp{z). (24) 

Again, it is straightforward to check that the diagonal component of this equation reduces 



to Eq.(|T2D. 



A. Case I: < 1 



Substituting Eq. (^) into (|T9D, we get, to leading order in 



Pab = 6ab-2j2 UaaUbaUbpUap siu^ ( ^Am^^ ) + 0{i 



a,f3 



L 



(25) 



Therefore, as expected, to this order, we obtain the standard expression for the transition 
probability well known for the three neutrino case. It is clear that in this scenario, solar 
and atmospheric data can be explained as in the usual three fiavour neutrino models by 
adjusting the spacing of the different ^q,'s. For the solar neutrino puzzle, one may either use 
MSW or VO solution depending on how much fine tuning one is willing to tolerate. 

For intermediate values where L/AE <C the leading corrections in ^ become 

important. Expanding Pab up to order by introducing Eq. 
standard oscillatory term we found 



1]) and neglecting the 
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where we have denoted 



Hz) 



— + c'iz) + s\z)-—c{z) 



(26) 



(27) 



Notice that the first term between parenthesis on Eq. (|26D contains the lower order correction 
to the standard survival probability in Eq. (pS]), while the second term, of oder will be 
only relevant for fiavour transitions. Taking the average on the last equations, and using 
that c = and + = 7r^/90, we get 



ab 



6. 



ab 
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TT 
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U aaUba^a 



(28) 



This last expression generalizes that presented in Ref. |Tl|, where only the contribution of 
^3 was assumed. 

It is clear that understanding LSND results in this case would require that first undergo 
a transition to the lower KK modes of the bulk neutrinos and then back to the Ve- One 
might hope that if we adjusted the extra dimension radius to be small enough ~ 0.2 — 2 



10 



eV, then one would get the right mass difference to fit LSND data. The key question then 
is to see whether the transition rate comes out right. For this purpose, we consider up to 
the lowest non- vanishing order in for P^e given in (pB]), which after using the unitarity of 
U and the fact that ^| — = 2Am^j^i?^ turns out to have the form 

P^e ^ sin^ 29,, {Aml^R') ' I{z) (29) 

with sin^26'^e = 4(?7^3[/e3)^- It is straightforward to check that P,e{.L = 0) = by using 
the identities c(0) = 7r^/6 and s(0) = 0. From the limits obtained by the CHOOZ and 
Palo Verde collaboration, we know that |f/e3p < 0.03 [|l^, and assuming |f/^3p = 0.5 for 
maximal mixing, the largest optimistic value for the mixing angle we may take is about 
sin^ 26'^e = 0.06. By fixing L/E as for LSND, I[z) becomes only a function of R (since 
z = L/2ER^). In figure 7 we have plotted I{R) versus R. From this figure we see that for 
a reasonable large R (~ 10 eV~^), the function I{R) has very small values already. The 
combined effect with the factor Arri^^^R^ < 1 will reduce P^e even more. Clearly, P^e will 
be maximal for the larger possible value of R. However, the largest allowed radius that 
permit us to still be confident in our approach is about 1/P ~ y^Am^^ ~ 0.06 eV. A 
numerical calculation with these inputs gives for LSND P^e = 3 x 10~^, which is one order 
of magnitude smaller than the observed anomaly in the u beam of LSND. Since higher order 
corrections in C,a to the probability are unlikely to introduce enough enhancement (a factor 
of 10 is needed), we conclude that this scenario yields a too small probability for z/^ — z/g 
transition to explain the LSND observations. 

B. Case II: ^1,2 < 1 < 6 

To proceed with this case, it is convenient to write the survival probability Paa in the 
following form: 

Paa= E (UaaUa^f p^^ + 2U!,J2 ULPa3 + U^3P33. (^O) 

a,l3=l,2 a=l,2 

The first term in the above equation can be written to leading order in the small ^'s as 
follows: 

E (t^aaMVa/3~ (l-f/a3)'-Sm'2^-Sin2(AA^2 \ +0(a (31) 
a,/3=l,2 ^'^^ ^ 

where sin^ 2^aa = '^^a\Ua2)^- This term by itself can explain the solar neutrino deficit. 
Clearly, the standard two neutrino oscillation expression is recovered to this order if we 
set f/e3 ~ to satisfy the bounds imposed by the reactor data [1^. This will make the 
contributions of the last two terms in the survival probability ( ^OD negligible and the first 
oscillatory term can then be used to solve the solar neutrino puzzle. One can of course use 
the MSW mechanism to solve the solar neutrino problem or the vacuum oscillation. The 
constraint on the parameter space is that the Amf2 appropriately adjusted. This does 
not impose any condition on the bulk radius and can be satisfied by the initial choice of 
parameters (the Yukawa couplings) in the theory. For instance, M* ~ 10 TeV can lead to 
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the MSW-type mass differences. Let us note that, if Am^2 ~ 10 ^ eV^, then, we will have 
|e?-e||~10-^i?^ (eV cmf. 

Now let us consider P^^. From Eq (0), we see that there are several contributions to the 
atmospheric neutrino deficit. First, there is the contribution of the towers labeled by ^i_2, 
which is oscillatory, although it can not be identified as in the usual u^^ —>■ Vr oscillations. 
Then, there is also the contribution induced by the term ^33 which is of the form of (|T2|). 
Finally, there is also a mixed term, p^s- 

To proceed with the full discussion, let us consider two cases: 
Case (i): ?7^3 = 

If f/^3 = 0, the last two contributions are removed and we get the survival probability 

P^^ ^ 1 - 4 {2U%Aml,R' + il) - c{z)^ (32) 

One might then hope that the oscillations into the lower KK modes with mass differences 
of about 1/R will do the job provided we choose 1/P ~ \J Am^j^ to match the data. In 
this case, the atmospheric muon neutrinos oscillate into the sterile neutrinos, a possibility 
which has its characteristic tests^. For this solution to work, one needs to assume ^1 ~ 1 
so that one gets maximal mixing. This in turn means that, to explain both solar and 
atmospheric neutrino data, an almost degeneracy ^1 ~ .^2 must be maintained. This alters 
our explanation of the solar neutrino deficit since, now, these new contribution to it (see 
the 0(^^) terms in Eq. (pH])), become more and more important and in fact of order one, 
making it hard to understand the solar neutrino deficit, since c{z) ~ 1. We will, therefore, 
consider this case as an unfavorable one for understanding the neutrino puzzles. 
Case (a): f/^3 7^ 

Turning to the case where f/^3 7^ 0, if we keep ^1^2 ^ 1 (to maintain our understanding 
of the solar neutrino data) then, the corrections of order C,i2 the transition probability 
become almost negligible and the dominant contribution to Pab then must come from ^3 
corrections. This yields 

Pab ~ {Sab - UazUbzf ' ^Ua^^ siu^ (^Am^,) 

+ 2Ua3Ub3 {Sab " UasUfi) C3{z) + ( f/a3 ^^63 ) ^3 • (33) 

Notice that P33 in the last equation has the same form as Eq. ([T^) , and the function C3 is 
the one defined above with the same sharp behaviour as P33. Then specializing Eq. (^) to 
the z/^ in the atmospheric case, we get 

P,, ^ (1 - Ul^y + 2t/i(l - UI3) C3{Z) + f/,^3 P33- (34) 

From our naive analysis in the previous section we may expect that this equation can account 
for the atmospheric data without too much trouble as long as ^|/P ~ 10~^ eV or so to 



^We realize that from an experimental point of view, this looks less hkely to be reahzed in 
nature [l^; we take a somewhat hberal view of the situation and still contemplate this as a viable 
possibility. 
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make the width of the slope larger than the experimental parameters and to avoid the over 



washing of the flux. Indeed, it has been checked numerically in reference |TT[] that the 
atmospheric neutrino data can be fitted in this case if U"^^ ~ 0.4 and ~ 0.02 eV. It is 
worth mentioning that a mixed explanation could be also possible, where the three towers 
contribute equally to provide atmospheric oscillations, but we will not discuss this case here. 

An important point to note however is that due to the features of the function 03(2), the 
atmospheric neutrino data will not exhibit the oscillatory behaviour that one would expect 
in the conventional two neutrino models. 

Lets turn now to LSND results. By taking that Ues ~ as suggested by solar neutrino 



and reactor data |]I5[, the transition probabihty reduces to P^e = I]q,/3=i,2 U^aUeaU^/sUe/s Pa/s- 
This reflects the fact that the same argument that suppresses the contribution of the third 
KK tower to Pee also does the same for P^e- As Ue3 = remove the contributions of the 
third tower, we may expand P^e to the lowest order by the same expression (|26| ) used in the 
previous case, which is now given as 

P^e-sm^2e^e {AmliR^fliz) (35) 

where now sin^ 2^^e = 4(f/^2f^e2)^; and I{z) as given in Eq. (^^. This resembles our former 
expression in (^). In order to estimate the magnitude of this contribution, note that the 
atmospheric neutrino fitting requires ^|/P ~ 10~^ eV. For ^| ~ 1 — 10, this implies that 
~ 10^ — 10^ eV~^. Thus if Aml^i corresponds to the MSW solution (large or small angle), 
then Aml^iR^ ~ 1 and one obtains P^e = sin^ 26'^e-^(2;)- From Fig. 7, we see again that for 
relevant values of R and L/E, the function I{z) takes very small values (~ 10~^), making 
the Pefj, very small. 

Notice that this argument is independent of the way we get the atmospheric deficit. We 
could also imagine, for instance, that a small value of Ues is allowed, and then calculate 
the leading correction to the above expression. However, it turns out to be of the form 
(U^sUes) J2a=i,2 U^aUea ^ai'^i^) ~ t)' "^^lere last term between parenthesis is already smaller 
than 10 ~^ by itself. 



C. Case III: 6 < 1 < ^2,3 
In this case, the transition probability can be written as 

Pab = {UalUblf Pll + 2UalUbl UaJJba Pal + UaJJbaUbpUap Pa/S- (36) 

a=2,3 ",/3=2,3 

The simplest possibility is to let the first term in above equation is be responsible for solar 



neutrino oscillations into bulk neutrinos as explained by Eq. (^ [10|. This requires that the 
radius be fixed to be about 1/P ~ 10~^ eV. In this case, to keep the z/g from mixing too 
much with the other neutrinos and generate further reduction of the survival probability for 
the solar neutrino, we choose Uei ~ 1. This essentially decouples the first tower from the 
others. As a simple approximation, if we assume that f/e,2,3 = 0, then clearly this suppresses 
the oscillations from z/^ into z/g, making it difficult to understand the LSND results. 

Moreover, in this scenario, even the explanation for atmospheric data seems to run into 
some trouble. As Uei ~ 1, it is not unreasonable to conclude based on orthogonality that 
U^i ~ 0. This implies that 
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-^MM — E i^fiaUf^/sY Pa/3, 
a,/3=2,3 



(37) 



where all pap are given as in Eq. (P^. As a rough approximation, if we assume that 
the partial transition probabilities are all almost of the same order (as they seem to be 
numerically), say P33, and use orthogonality once more to get P^^ ~ P33. Therefore, we 
get maximal contribution from the large deficit generated by P33. A hierarchical pajs will 
not help with this, since either we get equal mixing among them or one of them has a 
dominant contribution. In any case, the leading order will combine both things, a large and 
fast developing slope, and a large mixing angle. Nevertheless, based on the results of [jlT|, 
where they found some scenarios where those two ingredients come together (although for 
small (^), it is still hard to rule out the scenario without a careful numerical analysis of the 
data. In any case, the simplest requirement to get a deficit not too large compared with the 
experimental data, imposes a tuning of the main parameters ^23/-^ ~ 10~^, which combined 
with the condition for R to understand the solar neutrino data fixes C23 ~ 10- 

In order to understand the LSND results, we must allow for small Ue2 and/or [/^i. Since, 
in the present scenario, the solar neutrino deficit is being explained purely by oscillation 
into the bulk, so, we get the following picture: The first tower contributions are too small 
for LSND, they are even smaller than the size of those considered on the previous cases, 
while the oscillations involving the other towers are mainly destructive due to large values 
of ^ involved. The conversion will not occur but for a small contribution related with the 
small part of those towers contained in z/g. However, as one may suspect, the contributions 
are all proportional to the mixing angles {Ufj^aUeaY, which can not be larger than 10~^. So, 
the question is whether the behaviour of pap can help. To analyze this point we notice that 

Pf^e^ E U^aUeaUe(3U^f3 (l-Pafi)- (38) 

a,/3=2,3 

Note that one can estimate the value of (1 — Pa/s) from the the information on atmospheric 
neutrino deficit as follows. First point is that all Pa/s's are of same order and therefore, 
one can write the P^e ~ {U^iUei)^{l — Psurv), where PsurviC) is the same function which 
for ( = Catmos gives the survival probability of muon neutrinos in the atmospheric data. 
Thus, PsurviCatm) ~ 0.5. The value of ( corresponding to the LSND case is however much 
smaller due to smaller oscillation distance; therefore to estimate P^e we must use the value 
of ( which is much smaller and find the corresponding Psurv A numerical evaluation gives 
(1 — Psurv) ~ 10~^. The next question is how large the mixing parameters are. In order 
not to make Pee in atmospheric neutrino data different from one, niust be much smaller 
than 1. On the other hand to fit LSND observations, we need to have U^i ~ 0.5. It therefore 
appears difficult to accommodate the LSND observations in this case. 

IV. CONCLUDING REMARKS 

The analysis of the present paper shows that in minimal models for neutrino masses in 
theories with large extra dimensions, it is not possible to get a simultaneous explanation of 
gross overall oscillations needed to understand the solar, atmospheric and the LSND data. 
Fitting the overall deficit in the atmospheric and solar neutrino data for various ranges of 
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the input parameters, the largest conversion probabihty, P^e, that we find, is around 10~^. 
This is too low to explain the LSND observations. One must therefore invoke new physics 
beyond the standard model in the brane P,p!7| or new physics outside the brane |1I4| for a 



simultaneous understanding of all observed neutrino data. 

The oscillation pattern is governed by the dimensionless parameters ^a- Three of them 
arise in the models under consideration and we get the following picture: 

(i) ^,2.3 <^ 1: in this case, solar and atmospheric neutrino data are understood as in the 
case of four dimensional models and the hope that the presence of the bulk neutrinos 
with an appropriate KK excitation scale could provide an understanding of the LSND 
data is not realized. 

(ii) ^1^2 ^ 1 < 'Cs- solar neutrino data is provided as in the four dimensional models but 
atmospheric data is explained by to Ubuik oscillation once the appropriate values of 



the bulk radius is chosen i.e. either ^ ^2 and l/i? ~ yArriatm or ^|/-R ~ 10~^. 
In this case, the atmospheric neutrino fiux does not oscillate as a function of L/E, a 
feature distinct from the conventional two neutrino oscillation picture. 

(iii) ^1 <^ 1 < ^2,3- Both, solar and atmospheric data are explained by z/ — > Ubuik oscilla- 
tions. Therefore 1/ E? ^ ^^Li ~ 10~^ eV with matter effects for solar and ~ 10. 
Again, the L/E behaviour of the atmospheric flux is not oscillatory. 

(iv) ^1,2,3 ^ 1- There is no explanation for solar neutrino data in this case. This parameter 
range is therefore ruled out. 
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FIG. 1. Eigenvalues as calculated from Eq. (Q) for various values of ^. For comparison with 
the low ^ limit discussed in the main text we show the deviation from the integer numbers. This 
picture depict how for large ^ the eigenvalues get shifted close to semi-integer numbers. Notice 
also that for small ^; Aq is clearly the only non integer eigenvalue. 
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FIG. 2. Mixing factors for the expansion (^) for different values of the parameter ^. To 
make the effect for large ^ visible we have introduced an arbitrary normalization where Nq = 1. 
This amplifies the value of l/iV„ by a constant scaling. Of course, the real values satisfy the 
normalization condition. This figure shows how, for large ^, the effective number of KK mass 
eigenmodes contributing into the weak eigenstate vl gets larger and those modes become almost 
equally suppressed, while for small ^ the lightest mode (n = 0) become the main component, with 
the decoupling of all other modes. 
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FIG. 3. Exact survival probability for versus L/2E in units of B? (as it is explicit in 
the argument) for three different values of ^. Dotted lines represent the continuos approximation 
discused in the main text. Note that only the low ^ limit has a periodic behaviour. 
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FIG. 4. Here we show an amplification of the survival probability for the = 0.1 case showed 
in figure 3. Note the large number of wiggles produced by the oscillation of consecutive levels in 
Eq. (pi). We also depict the continuos limit (dotted line) for comparison. 
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FIG. 5. In this figure we plot the functions p{C) s-^icl i?(C) that describe the slope of the survival 
probability Pgurv (which we also show) in the continuos limit where ^ > 1. The argument is the 
dimensionless parameter C = {tt^"^ /2R)yJ L/2E defined in the main text. 
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FIG. 6. Dependance of the survival probability as a function of ^ for the characteristic 
experimental values of L/AE, with R = 0.2 mm. 
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FIG. 7. The function I{R) versus R for the experimental ratio L/2E as in LSND. The window 
on the upper right shows an amplification of the region for large R, where I{R) < 10"'^. 
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